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Abstract 



^P I Quantum field theories require a cutoff to regulate divergences that result from local 

interactions, and yet physical results can not depend on the value of this cutoff. The renor- 
^S . malization group employs a transformation that changes the cutoff to isolate hamiltonians 

^P I that produce cutoff- independent eigenvalues. The similarity renormalization group is based 

on similarity transformations that regulate off-diagonal matrix elements, forcing the hamil- 
*pH, tonian towards a band-diagonal form as the cutoff is lowered. This avoids pathologies that 

0^\ plagued tradition transformations acting on hamiltonians, making it possible to produce 

a well-behaved perturbative approximation of renormalized hamiltonians in asymptotically 
free theories. We employ a simple two-dimensional delta function example to illustrate this 
new renormalization technique. 






X 



1 Introduction 

Early attempts to combine quantum mechanics and special relativity led to the consideration 
of local interactions, which are consistent with causality and avoid signals that propagate 
faster than light. Local interactions lead to divergences in perturbation theory, whose dis- 
covery caused some of the best theorists in the world to question the foundations of quantum 
mechanics. Eventual successes at fitting precise atomic experimental data led to the univer- 
sal acceptance of renormalization recipes that were acknowledged to make little sense Q. 
Initially the perturbative renormalization of QED required theorists to match perturbative 
expansions in powers of a bare and physical electronic charge [^ , but the bare charge clearly 
diverges logarithmically in QED and the success of an expansion in powers of such a coupling 
was mysterious at best [^. 

The first steps towards making sense of renormalization theory were taken in the 1950's 
with the invention of the perturbative renormalization group |^ ^, ^, 0, P, although serious 
investigators found the theory was still plagued by non-convergent sums because QED is not 



asymptotically free. The development of Wilson's renormalization group formalism y, 10 1 
and the discovery of asymptotic freedom [|l^] allowed physicists to produce a logically rea- 
sonable picture of renormalization in which perturbative expansions at any high energy scale 
can be matched with one another, with no necessity to deal with intermediate expansions in 
powers of a large parameter. 

In this pedagogical article we take advantage of the fact that the divergences in field 
theory result entirely from local interactions. To understand the most important aspects 
of renormalization theory requires only a background in nonrelativistic quantum mechanics, 
because as has been long known the divergences of field theory are directly encountered when 
one tries to impose locality on the Schrodinger equation. In this case the interactions we 
consider that are local at all scales are delta functions and derivatives of delta functions. 
These divergences can be regulated by the introduction of a cutoff, and the artificial ef- 
fects of this cutoff must be removed by renormalization. The simplicity of the one-body 
Schrodinger equation makes it possible to renormalize the theory exactly, disentangling the 
effects of locality from the complicated many-body effects and symmetry constraints encoun- 
tered in realistic field theories. There is a large literature on the subject Il^]-p9[, primarily 
pedagogical. 

The similarity renormalization group (SRG) is a very recent development invented by 



Stan Glazek and Ken Wilson [^, 31], and independently by Franz Wegner [^. We do not 



review the applications of this method, which are growing in number. 

In the SRG, as in Wilson's original renormalization group formalism |33, 34], transfor- 
mations that explicitly run the cutoff are developed. These transformations are the group 
elements that give the renormalization group its name. 

In his earliest work ]^, ^J] Wilson exploited a transformation originally invented by 
Claude Bloch [p3]. It uses a cutoff on the states themselves, and as the cutoff is lowered, 
states are removed from the Hilbert space. If the hamiltonian is viewed as a matrix, these 
cutoffs can be seen as limiting the size of this matrix and the transformation reduces this 
size, as illustrated in Fig. la. Wilson introduced a rescaling operation to allow transformed 
hamiltonians to be compared with initial hamiltonians, despite the fact that they act on 
different spaces; however, the Bloch transformation is ill-defined and even in perturbation 
theory it leads to artificial divergences. These divergences come from the small energy dif- 



/ 



\ 







I 







/ 




(a) 



(b) 



Figure 1: Two ways to run a cutoff on free energy. In (a) a cutoff on the magnitude of the energy 
is lowered from the sohd to the dashed hues, with problems resulting from the removed shaded 
region. In (b) a cutoff on how far off diagonal matrix elements appear is lowered from the dashed 
to the solid lines. 



ferences between states retained and states removed by the transformation, and they appear 
in the form of small energy denominators in the perturbative expansion of the transformed 
hamiltonian. These small energy denominator problems led Wilson to abandon the hamilto- 
nian formulation of field theory in favor of path integral formulations, but the virtues of the 
hamiltonian formulation over the path integral formulation for many problems remains. 

The breakthrough provided by the SRG is that the transformations are typically unitary, 
making them well-defined, and they run a cutoff on energy differences rather than on indi- 
vidual states, as illustrated in Fig. lb. Again viewing the hamiltonian as a large matrix, 
these cutoffs limit the off-diagonal matrix elements and as they are reduced the hamiltonian 
is forced towards diagonal form. The perturbative expansion for transformed hamiltonians 
contains no small energy denominators, so the expansion breaks down only when interactions 
become sufficiently strong, in which case perturbation theory should fail in any case. 

Although the SRG has not yet been applied to a wide range of problems, it may be an 
important new tool both for attacking field theories and non-relativistic many-body problems. 

When the SRG is used with coupling coherence |36, 37 1, which we explain below, it 



allows us to construct effective theories with the same number of free parameters as the 
underlying 'fundamental' theory. For the delta-function example there is one fundamental 
parameter, the strength of the regulated delta- function as the cutoff is removed. In the SRG 
with coupling coherence, there is only one fundamental coupling and all new couplings are 
perturbative functions of the fundamental coupling that are given by coupling coherence. 



It is the renormalization group flow of tlie added couplings, and a boundary condition that 
they vanish when the fundamental coupling is taken to zero, that fixes their dependence on 
the fundamental coupling. 

The examples we use in this article do not illustrate non-Gaussian fixed points, so their 
scaling properties are driven by naive dimensional analysis. However, we will see that even in 
these cases scaling behavior of effective hamiltonians derived using a perturbative similarity 
renormalization group can be very complicated. We will see that in the perturbative SRG 
there are errors arising from the approximate treatment of the fundamental running coupling 
and the approximate treatment of the relation between this coupling and the new couplings 
of irrelevant operators. 

In a realistic calculation the marginal coupling, which corresponds to the strength of the 
regulated delta function, would be fit to data. In order to clearly illustrate the logarithmic 
errors that result from using the perturbative SRG equations, we approximate this marginal 
coupling in this article rather than renormalizing it nonperturbatively by fitting data. The 
strengths of the irrelevant operators, which correspond to derivatives of the regulated delta 
function, are approximated using expansions in powers of the approximate running coupling 
that are fixed by coupling coherence. The approximate running coupling differs from the 
exact running coupling by inverse powers of logarithms of the cutoff, and the error analysis 
for the binding energy displays the resultant inverse logarithmic errors in addition to power- 
law errors seen in all approximate renormalization group calculations. In addition there are 
errors in the strengths of the irrelevant operators resulting from using a truncated expansion 
in powers of the running coupling and an approximate running coupling, both of which 
introduce inverse logarithmic errors in addition to the power-law errors normally seen. 

The utility of the renormalization group rests on our ability to accurately determine 
and control the magnitude of errors resulting from the artificial cutoff. For perturbative 
calculations this issue is not critical, but in all field theories (and in our example) a scale is 
reached where the coupling becomes large and a non-perturbative calculation must be done. 
The renormalization group allows us to eliminate as much perturbative physics as possible 
(i.e., lower the cutoff as far as possible in an asymptotically free theory), so that the essential 
degrees of freedom that couple non-perturbatively can be isolated. 

2 Similarity Renormalization Group 

In this section we review the general formulation of the SRG developed by Glazek and Wilson 



|30| , 31] and a specific transformation developed by Wegner ^^. The reader may wish to 



skip the general formulation on a first reading. 

2.1 Glazek- Wilson Formulation 

Consider a system described by a hamiltonian written in the form 

H = h + V , (1) 

where h is the free hamiltonian and V is an interaction. 

In general, the hamiltonian can couple states of all energy scales and such couplings can be 
a source of ultraviolet divergences. The goal of the SRG is to obtain an effective hamiltonian 



in which the couphngs between high and low-energy states are removed, while avoiding 
any problems from small energy denominators in effective interactions. The procedure is 
implemented by a unitary transformation that generates effective interactions that reproduce 
the effects of the eliminated couplings. The effective hamiltonian cannot produce ultraviolet 
divergences at any order in perturbation theory as long as its matrix elements are finite. 
In our discussion we will use the basis of eigenstates of the free hamiltonian, 

h\i >= ei\i > . (2) 

We start by defining a bare hamiltonian, H\, regulated by a very large cutoff A (here with 
dimensions of energy) on the change in free energy at the interaction vertices, 

Ha = h + VA, (3) 

Va = /aFa, (4) 

Va = v + Hf, (5) 

where /a is a "similarity function", Va is defined as the reduced interaction and W^ are 
counterterms that must be determined through the process of renormalization in order to 
remove A dependence in physical quantities. 

The similarity function /a regulates the hamiltonian by suppressing matrix elements 
between free states with significantly large energy difference and acts in the following way: 

<^I/a^aIj> = ei'Jii + /A(ei-ej) <^|Fa|J > 

= ^i^ij+ fMi^Mj- (6) 

Typically, the similarity function is chosen to be a smooth function satisfying 

(0/A(ei - ^i) -^ 1' when je, - ej\ « A , 

{ii)f^{ei - ej) -^ 0, when |ej - ej\ » A . (7) 

In several papers the similarity function has been chosen to be a step function. Although 
useful for doing analytic calculations, such a choice can lead to pathologies. 

The similarity transformation is defined to act on the bare regulated Hamiltonian, Ha, 
lowering the cutoff down to a scale A: 

Hx = U{X,A)HaU\X,A). (8) 

The renormalized Hamiltonian can be written in the general form 

Hx = h + Vx, (9) 

Vx = /a Fa. (10) 



The transformation is unitary, so Ha and Hx produce the same spectra for observables. Also, 
if an exact transformation is implemented, the physical predictions using the renormalized 

TCt 

'a 
independent of A as A ^ oo. 

The unitarity condition is given by: 



Hamiltonian must be independent of the cutoff A and H'^ is chosen so that they also become 



U{X,A)U\X,A) = U\X,A)U{X,A) = 1 . (11) 



The similarity transformation U can be defined in terms of an anti-hermitian operator Tx 
4 



(T| = —Tx) which generates infinitesimal changes of the cutoff energy scale, 



C/(A,A) = Texpl f Ty dX'j , 



(12) 

where T orders operators from left to right in order of increasing energy scale A'. Using 

Tx = U{X, A) ^^^^ = yj^ f/t(A, A) , (13) 

dX dX 



and the unitarity condition Eq. (11), we can write Eq. (^) in a differential form, 

dH^ 



dX 



[Hx,Tx] . (14) 



This is a first-order differential equation, which is solved with the boundary condition 
Hx\x->A — -^A- The bare Hamiltonian is typically given by the canonical Hamiltonian plus 
counterterms that must be uniquely fixed to complete the renormalization. 

The operator Tx is defined by specifying how Vx and h depend on the cutoff scale A. For 
simplicity in this article, we demand that h is independent of A, although this may not lead 
to an increasingly diagonal effective hamiltonian in all cases. We also demand that no small 
energy denominators can appear in the hamiltonian. These constraints are implemented by 
the conditions 

^ - IV^-T.]. (16) 

To obtain the renormalized Hamiltonian perturbatively, expand 

— —(1) —(2) , V 

Vx = Vx +Vx +■■■ , (17) 

rA = rf +rf + ••• , (is) 

< = <" + <^' + --- (19) 

where the superscripts denote the order in the original interaction, V. A general form of 
these effective interactions is 



r-i °° /■A 

Vx = - E '^O+M / dX' [Fif ,Tir] +<'^* , (20) 

^ 1 — 1 -^A 



_(i) 

for i = 2, 3, • • •, with V x = v. For instance, the explicit form of the second-order effective 

_(2) 

interaction V x is 

^S=E^^^J^+¥^ +<.^ (21) 



A,;. A 



■jk 



where 

(AA) _ [^ ,w „ dfx'ki /r,o\ 

^'''^ Jx '^^ ~d>r ' ^ ' 

^ij = Cj - Gj . (23) 

(n) ct 

The counterterms Hj^ ' can be determined order-by-order using the idea of couphng- 
coherence |36, 5^. This is implemented by requiring the hamiltonian to reproduce itself 



in form under the similarity transformation, the only change being explicit dependence on 
the running cutoff in the operators and the implicit cutoff dependence in a finite number 
of independent running couplings. All other couplings depend on the cutoff only through 
their dependence on the independent couplings. In general, we also demand the dependent 
couplings to vanish when the independent couplings are taken to zero; i.e, the interactions are 
turned off. If the only independent coupling in the theory is ax, the renormalized hamiltonian 
can be written as an expansion in powers of this coupling: 

Hx = h + axO^^^ + ai0(2) + . . . . (24) 

In this way, the effective hamiltonian obtained using the similarity transformation is com- 
pletely determined by the underlying theory. The procedure can be extended to arbitrarily 
high orders, although it becomes increasingly complex both analytically and numerically. 

2.2 Wegner Formulation 

The Wegner formulation of the SRG [p2| is defined in a very elegant way in terms of a flow 



equation analogous to the SRG Equation in the Glazek- Wilson formalism |3y, 31], 



^ = [H.,n. (25) 

Here the hamiltonian Hs = h + Vg evolves with a flow parameter s that ranges from to oo. 
The flow-parameter has dimensions 1/ (energy)^ and is given in terms of the similarity cutoff 
A by s = 1/A2. 

In Wegner's scheme the similarity transformation in defined by an explicit form for the 
generator of the similarity transformation, Tg = [Hs,Hq\, which corresponds to the choice 
of a gaussian similarity function with uniform width. In the original formulation, Wegner 
advocates the inclusion of the full diagonal part of the hamiltonian at scale s in Hq. For a 
perturbative calculation of Hg, we can use the free hamiltonian, Hq = h. With this choice, 
the flow equation for the hamiltonian is given by 

^ = [Hg,[HgM]- (26) 

The reduced interaction, V gij (the interaction with the gaussian similarity function fac- 
tored out) is defined by 

Vsij = A,, Vgi^ , (27) 

/., = e-^^- . (28) 



Assuming that the free hamiltonian is independent of s, we obtain the flow equation for the 
reduced interaction, 

^ = ^ (A,fc + A,fc) Vs^k Vsk^j e-2^^-^- , (29) 

^* k 

where we use A|- — A^^ — A^^ = —2/S.ik/S.j].. We should emphasize that this is an exact 
equation. 

To solve this equation we impose a boundary condition, Hs\^^^ = ^sq- Then, we make 
a perturbative expansion, 

Vs = V, +V, +••• , (30) 

where the superscript implies the order in the bare interaction V s ■ It is important to 
observe that counterterms are implicit in the bare interaction and can be determined in the 
renormalization process using coupling coherence. 
At first order we have 

—r^ = , (31) 

as 
which implies 

vZ = Vs,^J, (32) 

where s is the final scale. Because of the dimensions of the flow parameter we have s > Sy, 
corresponding to a smaller cutoff. The "no cutoff limit" corresponds to s^ — > 0. 
At second order we have 

_(2) 

-^ = Yl (A^fe + A,fe) Vs^^k Vs^^kj 6-2^^^'=^'= . (33) 

k 

Integrating, we obtain 



—(2) _1„_ _ / I 1 \ 



-2soA,fcAjfe _ -2sAikAjk 



(34) 



By construction, the Wegner transformation is unitary and avoids small energy denom- 
inators. The Wegner transformation is one of the Glazek- Wilson transformations, with the 
similarity function chosen to be f^^. = e~ ^r . 

2.3 Strategy 

In our applications of the SRG we use Wegner's transformation. The renormalized hamilto- 
nian for the non-relativistic delta-function potential in D-dimensions is given by 

H,{p,p')=p'6^^\p-p') + e~^'^lf^ fF«(p,p') + yf (p,p') + ...l , (35) 



where 

^i'^(P'P') = -^^ (36) 

yf(p,p') = 4,Ff)(p,p'), (37) 

yi")(p,p') = al,F(^\p,p'). (38) 

Here A is a momentum cutoff (as opposed to the energy cutoff discussed above) related to 
the flow parameter by s = 1/A^. The index i denotes the order of the calculation for the 
running coupling. 

The renormalized hamiltonian can be used to compute eigenvalues and eigenstates. Since 
the hamiltonian is derived perturbatively we expect cutoff dependent errors in the observ- 
ables. Formally, we can regroup the terms in the renormalized hamiltonian and write it as a 
momentum expansion, and the expansion parameters are analytic functions of the running 
coupling a\. Expanding the operators -Fs (p,p') in powers of p^/A^ we obtain 

Fs '{P,P)=Z0 + Z2 ^ + Z4 ^^4 + ^4 -^^ + ■■■ , (39) 



(P^-P'^)^ 



where the Zi s are constants. Regrouping the terms we obtain 

Hxip,p')=pH^''\p-p')+e 
where 



(p^ +p'^) 
go{ax) + g2{ax) — ^ \- ■■■ 



(40) 



gi{ax) = aiax + bial + ... . (41) 

We can identify three interdependent sources of errors in the perturbative similarity 



renormalization group when the hamiltonian given by Eq. (40) is truncated and used to 
compute a physical quantity: 

a) errors introduced by the truncation of the hamiltonian at a given order in p^ /}?; 

b) errors introduced by the truncation of the hamiltonian at a given order in the running 
coupling aA,i, which correspond to the use of an approximation for the functions gi] 

c) errors introduced by the approximation for the running coupling ax,i. 



In the actual calculation using the hamiltonian given by Eq. ( |35| ) errors of type (a) do not 
appear directly because we do not truncate the operators that appear in the hamiltonian. 
However, errors of type (b) can be understood as coming from approximating the couplings 
in front of the operators in Eq. (HO). Errors of type (c) appear in our calculations only 
because we do not fit the canonical coupling to data at each scale, but fix it at a given scale 
and evolve it perturbatively from that scale. The strategy we would use for a realistic theory 
{e.g., QED and QCD) is the following: 

1) Obtain the renormalized hamiltonian using the similarity transformation and coupling- 
coherence, truncating the hamiltonian at a given order in powers of Oix,i- 

9 



2) Fix the coupling ax by fitting an observable {e.g., a bound-state energy). 

3) Evaluate other observables {e.g., scattering phase shifts). 

As pointed out before, the evaluation of scattering observables with the similarity hamil- 
tonian with standard techniques is complicated and so in our examples we focus on the bound 
state errors. We fix the coupling at some scale using a given renormalization prescription and 
use the flow-equation to obtain the coupling as a function of the cutoff A to a given order. 
We then perform a sequence of bound-state calculations with better approximations for the 
hamiltonian such that the errors in the bound-state energy are systematically reduced. Once 
the sources of errors are identified, it becomes relatively simple to analyze order-by-order 
how such errors scale with the cutoff A. In principle, to completely eliminate the errors 
proportional to some power m in the momentum expansion we should use the similarity 
hamiltonian with the exact running coupling (renormalized to all orders) and include the 
contributions up to 0{p"^/X"^) coming from all effective interactions (all orders in ax). Some 
details of this scaling analysis are presented later for the specific examples we work out. We 
emphasize again that in a realistic calculation we would fit the coupling ax to an observable. 
This nonperturbative renormalization eliminates the dominant source of errors we display 
in SRG calculations in this paper. We choose to renormalize the coupling perturbatively 
in this paper because the only observable we compute is the single bound state energy of a 
delta-function potential, and fitting this energy would prevent us from displaying errors. 

3 Two- Dimensional Delta-function Potential 

We now consider the case of two nonrelativistic particles in two dimensions interacting via 
an attractive Dirac delta-function potential. The Schrodinger equation for relative motion 
in position space (with ^ = 1), can be written as: 

- V^^(r) - Qo 5^^\r) *(r) = E 'f (r) . (42) 

Both the delta-function potential in two dimensions and the kinetic energy operator scale as 
1/r^, therefore, the coupling ao is dimensionless. As a consequence, the hamiltonian is scale 
invariant {i.e., there is no intrinsic energy scale) and we can anticipate the presence of loga- 
rithmic ultraviolet divergences, analogous to those appearing in QED and QCD. The problem 
requires renormalization. In this subsection we present the standard method that produces 



an exact solution analytically, using simple regularization and renormalization schemes [18|. 



3.1 Exact Solution 

We start with the Schrodinger equation in momentum space, 

/ Hp) -J^ J d'q ^(q) = E ^p) , (43) 

where $(p) is the Fourier transform of the position space wave- function, 

$(p) = ^ f d^r *(r) e-^P"" . (44) 

27r J 

10 



As a consequence of scale invariance, if there is any negative energy solution to Eq. (^) 
then it will admit solutions for any E < 0. This corresponds to a continuum of bound 
states with energies extending down to — oo, so the system is not bounded from below. By 
rearranging the terms in the Schrodinger equation we obtain 

Hp) = ^ ^^"^ (45) 

where ^(0) is the position space wave- function at the origin and iiJo > is the binding energy. 
To obtain the eigenvalue condition for the binding energy, we can integrate both sides of 
Eq. Q: 

' = ?. r ""'WTE^y W 

The integral on the r.h.s. diverges logarithmically, so the problem is ill-defined. 

The conventional way to deal with this problem is renormalization. First, we regulate 
the integral with a momentum cutoff, obtaining 



ao 



^ . 1 




so that 



1 = ^/ dpp , , ^ =^ln{l + -^] , (47) 

Eo = ^jg . (48) 

e "0 — 1 

Clearly, if the coupling oq is fixed then E'q ~^ oo as A — > oo. In order to eliminate the 
divergence and produce a finite, well-defined bound state we can renormalize the theory by 
demanding that the coupling runs with the cutoff A in such a way that the binding energy 
remains fixed as the cutoff is removed: 

«o ^ "A = —J- ^jy • (49) 

The dimensionless renormalized running coupling a^ that characterizes the strength of 
the interaction is therefore replaced by a new (dimensionful) parameter £"0 > 0, the binding 
energy of the system. This is a simple example of dimensional transmutation jSg]: even 
though the original "bare" hamiltonian is scale invariant, the renormalization procedure leads 
to a scale that characterizes the physical observables. Note that Eq can be chosen arbitrarily, 
fixing the energy scale of the underlying (renormalized) theory. It is also interesting to 
note that the renormalized running coupling a\ vanishes as A — > 00 and so the theory is 
asymptotically free. 

This renormalized hamiltonian can be used to compute other observables. The usual 
prescription for the calculations is to obtain the solutions with the cutoff in place and then 
take the limit as the momentum cutoff is removed to 00. If an exact calculation can be im- 
plemented, the final results should be independent of the regularization and renormalization 
schemes. As an example, we calculate the scattering wave function, 
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where k = \/E. Integrating both sides over p with a cutoff A in place, we obtain 

-1 



m = ^ 



47r \ — K^ — le 



thus, 



OA ^(0) 



27r 



1 / A2^ 

4tt \ Eq 



if A2 ^ 
— In 1 + — ^ 

47r \ — K^ — ze 



n -1 



In the limit A ^ co we obtain: 



ttA ^(0) 



In 



(51) 



(52) 



(53) 



I TT 



The resulting scattering wave function is then given by 



<^kip) = s^'Hp-^) + ,^ 



2tt [p^ — k"^ — ie) 



-f 



n -1 



I vr 



(54) 



It is important to note that only S-wave scattering occurs, corresponding to zero angular 
momentum states. For the higher waves the centrifugal barrier completely screens the delta- 
function potential and the non-zero angular momentum scattering states are free states. 

The same prescription can be used to evaluate the T-matrix or the K-matrix. For the 
T-matrix, the Lippmann-Schwinger equation with the renormalized potential is given by: 



r(p,p';A;)=y(p,p')+y' A 






(55) 



^2 _ g,Z _|_ ^g 

Since only S-wave scattering takes place we can integrate out the angular variable, obtaining 

fA 
/O 

where 



r('=°)(p,p'; k) = V(^='\p,p') + dqq ^V ''':".' T(^='Hq,p'; k) 



^2 _ g,2 _j_ ^g 



V^'='\p,p') 



OA 

' 2tt 



(56) 
(57) 



The Lippmann-Schwinger equation for the "on-shell" T-matrix is given by: 



^,.„,,„ _«,_«,,„.„,,,, ^ ,,,_^ 



+ ie 



(58) 



Solving this equation and taking the limit A ^ 00, we obtain the exact "on-shell" T-matrix: 

2 



To(fc) 



l-'f 



(59) 



I vr 



Here and in what follows we drop the superscript and use the subscript to denote the exact 
quantities. 

In the same way, the S-wave Lippmann-Schwinger equation for the K-matrix is given by 



K{p,p'-k) = V{p,p')+rJ^ dq q J^ K{q,p';k) , 



(60) 
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and the exact "on-shell" K-matrix is given by 



^o(/c) = -— ^. (61) 

The "on-shell" K-niatrix and T-niatrix are related by 

K^i^k) = ^°^^^ . (62) 

Using either 

2 fc 1 

k cot 5Q{k) -ik = — -— - , (63) 

vr lQ{k) 

or 

A; cot <5o(/c) = — — -, (64) 

vr Ao(k) 

we can obtain the exact phase-shifts: 

cot 5o = - In f 1^ 1 . (65) 

vr V^o/ 

3.2 Similarity Renormalization Group Approach 

In the two-dimensional case the canonical hamiltonian in momentum space with a delta- 
function potential can be written as 

i/(p,p') = /i(p,p') + np,p'), (66) 

where /i(p, p') = p^S^"^' (p— p') corresponds to the free hamiltonian and l^(p, p') = — ao/(2vr)^ 
corresponds to the Fourier transform of the delta-function potential. 

Integrating out the angular variable, the flow equation obtained with Wegner's transfor- 
mation in terms of matrix elements in the basis of free states is given by 

^Y^kl^l = -ip^-p'^f Vs{p,p') - r dkk (2A:2 -/ -p'2) V,{p,k) Vs{k,p') . (67) 
as Jo 

In principle, we can set the boundary condition at s = (no cutoff), i.e, 

Hs=o{p,p') = H{p,p') =p'5^'\p-p') - ^ . (68) 

However, the hamiltonian with no cutoff produces logarithmic divergences, requiring renor- 
malization. As we will see, the boundary condition must be imposed at some other point, 



leading to dimensional transmutation |38|. The reduced interaction Vs{p,p') is defined such 
that 

Vs{p,p') = e-<P"-P'"'>"Vs{p,p'). (69) 

Assuming that h is cutoff independent we obtain the flow equation for the reduced interaction, 

ds Jo 

xVs{p,k)Vs{k,p'). (70) 
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This equation is solved using a perturbative expansion, starting with 



'2lT 



(71) 



We assume a couphng-coherent solution in the form of an expansion in powers of as/2Tr, 
satisfying the constraint that the operators Fg {p,p') vanish when p = p' = 0, 



Vs{p,p') 



as 
'2tt 



E 

n=2 



as 

27r 



Fi")(p,pO 



(72) 



Note that the expansion parameter is as/27r. 

Using the solution Eq. (^) in Eq. (|70| ) we obtain 



ds 



1 dos ^ 1 



?{n), 



n 



ar'^-^Fi->(p,p')+a^'J^-Ml 



dk k (2fc2 _ p2 _ p'2) e-2.b'p'»+''-''(P^+P' 



S + Ef£) n'-'fe^ 



27r 



n=2 



m=2 



S + E S ^i™H^'^') 



(73) 



This equation is solved iteratively order-by-order in as/2TT. Again, if as/2TT is small the 
operator V} {p,p') can be identified as the dominant term in the expansion of Vs{p,p') in 
powers of p and p' . In the D = 2 case this operator corresponds to a marginal operator 
(since the coupling is dimensionless and there is no implicit mass scale). The higher-order 
terms correspond to irrelevant operators. 
At second-order we have 



1 da. 1 2 dFP{p,p') 



27r ds (27r)2 



ds 



-a. 



If\p,p'), 



(74) 



where 



If\p.p') 



(2^)2 Jo 



dk k {2k'' -p'-p'') e-2.[pV^+fc*-fc2(p2+p 



1 e"^'' P '''' 



(27r)2 4s 
The equation for a. is obtained by taking the limit {p.,p') -^ 0, 



where 



Integrating Eq. (ffq) from sq to s 



1 dos 

2tt ds 

If (0,0) 
as,2 



a. 



If) (0,0) 



1 1 



(27r)2 4s 



a 



so 



1 _ !i£CL In f ^ V 



(75) 

(76) 
(77) 

(78) 
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In terms of the cutoff A we obtain 



a\,2 



"Ao 



l + ^ln i 



(79) 



In principle, knowing the value of a^g for a given sq we can determine the running coupling 



Us for any s. Since we cannot choose sq = (Aq = oo), to use Eq. (75) we must specify a 
renormalization prescription that allows us to fix the coupling at some finite non-zero value 
of sn. We discuss this issue in detail later in this subsection. 
The equation for Fs {p,p') is given by 



1 dFp\p,p' 
(27r)2 ds 

Integrating from sq to s we obtain 



lPiO,0)-lP{p,p') 



(80) 



fP{p,p') 



ds' 



I _ g-2s' p^p'- 



so 



is' 



In 



So 



Ei(-2s p^ p'-^) + Ei(-2so ^ P ) 



.(2) 



Insisting that Fg {p,p') = when p,p' = we obtain 



fP{p,p') 

At third-order we have 



1 



7 + ln(2s p^ p'^) - Ei(-2s p^ p'^) 



(81) 



(82) 



1 das 



1 2 dFp\p,p') 



27r ds (27r) 



2 s 



ds 



+ 



da^ 



(27r)2"^ ds ' ^^'^^+(2^)3"^ 
-a2/(2)(p,p') + a^lf(P,l''), 



1 3 dFi''\p,p') 



ds 



(83) 



where 



iPip^p') 



(27r)3 7o 



dit it {2k^ -p'-p'^) ^-2s[p\'^+k^-kHp-+p'^)] 



X FP(p,A:)+FP(fc,p'y 
In the limit p,p' ^ we obtain: 



^^ = «^/?H0,0)-a3j(3)(0,0), 



where 



If) (0,0) 



(27r)3 Jo 



dk 2^^6-2^'=' 



FP(0,A:)+Fi2)(A;,0)' 



(84) 



(85) 



(86) 



.(2), 



n(2) 



Since Fg {k, 0) = F^ (0, k) = 0, the term proportional to af in Eq. ( |8g) is zero 



15 



The equation for Fg {p,p') is given by 



p(3). 



ji^'-^^ = -I ^^^'M Fi^\p,p')^ii^\p,p') . 



(87) 



n{3), 



To obtain Fg {p,p') the integrals in k and s must be evaluated numerically. 
At fourth-order we obtain 



1 das 

2-ir ds 



+ 



1 2 dFP{p,p' 



{2^f 



a^ 



a„ 



ds 
dFP{p,p') 



+ (2^"^^^^^'^^^'^') 



.(4), 



, 3 2 ^"s 77(3)/ /x , 1 4 '^^'^ (P,y) 



(27r)3 ^^ ds 



(88) 



where 



li'Hp,/) 



(27r)4 Jo 



'dA; fc {2^ -p'-p") ^-2slp^P''+k^-kHp'+p")] 



fP{p, k) + FP{k,p') + Fi2)(p, A;)Fi2)(A;,/^ 



In the limit p,p' — > we obtain: 

1 (ia^ 
27r ds 



a, 



/p)(o,0)-af/f)(0,0), 



where 



7(4) (0,0) 



1 

2^ 7o 



dA; A; 2A;2e-2^'=' 



Fi3)(0,fc)+Fi3)(fc,0) 



For dimensional reasons Eq. (|90| ) takes the form 

^2 



das 



_ 2 _ :^ 4 
s s 



(89) 



(90) 



(91) 



(92) 



where B2 = -^ and B4 can be obtained by evaluating ll (0,0) (numerically) for s = 1. In 
terms of the cutoff A we obtain: 



da\ 
~d\ 



1 



1 



^4 ^2 «! - -^ -B4 a\ . 



54 a^ 



(93) 



Integration of Eq. (pQ ) leads to a transcendental equation which is solved numerically in 
order to obtain the running coupling ag^i- 

Qualitatively, the errors are expected to be a combination of inverse powers of A and 
powers (or inverse powers) of ln(A) coming from the perturbative expansion in powers of a\ 
for the coefficients of the irrelevant operators and the perturbative approximation for a\. 

As pointed out above, to completely determine the renormalized hamiltonian at a given 
order we need to specify the coupling at some scale Aq- The simplest way is to choose a 
value for the 'exact' a^g- Formally, this fixes the underlying theory; i.e., if we had the exact 
hamiltonian (to all orders) we could obtain the exact values for all observables. However, since 
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the hamiltonian is derived perturbatively and must be truncated at some order in practical 
calculations, we can only obtain approximate cutoff-dependent results for the observables. 
Moreover, in this case the errors cannot be directly evaluated, since the exact values remain 
unknown. As an example, we calculate the bound-state energy choosing ax^ = 1.45 at 
Ao = 100. In Fig. 2 we show the binding-energy as a function of the cutoff A using the 
following approximations for the interaction: 



(a) marginal operator with coupling (oao), 

/ 2_ /2^2 

vx{p,p') = -'^^'^^ ; 

(b) marginal operator with running coupling renormalized to second-order (0^,2)1 

QA,2 _(£!V!)! 



Vx{p,p') 



27r 



(94) 



(95) 



(c) marginal operator plus second-order irrelevant operator with running coupling renormal- 
ized to second-order {a\2iF^ 



A ^ 



yx{p,p) 



a\,2 (a\^2 



Ff\p,p') 



P 



27r ■ V 27r 
(d) marginal operator with running coupling renormalized to fourth-order (aA,4) 



(96) 



y\{p,p) 



2tt 



— p — 



(97) 



(e) marginal operator plus second-order irrelevant operator with running coupling renormal- 

(2) 
ized to fourth-order a\^i,F^ ), 

Vx{p,p') = 



aA,4 fax^A 



2tt 



V 27r 



Fi'\p,p') 



, 2 /2n2 
P 



(98) 



We see that as the approximation is improved the cutoff dependence is reduced. As 
A ^ 00 all curves should approach the same binding-energy, which corresponds to the exact 
value, and as A becomes small the perturbative approximation breaks down. 

A similar prescription is to find Oaq at Aq that produces a given binding-energy, Eq. Since 
the fitting is implemented using a truncated hamiltonian, this oaq is an approximation that 
becomes more accurate if we use a larger Aq and/or include higher order operators. Although 
in this case we can evaluate the errors, the scaling analysis becomes complicated as A ^ Aq 
because at this point we force the energy to be the exact value and so the error is zero. 

As an example, we calculate the bound-state energy when the coupling is fixed at Aq = 100 
to give £"0 = 1. In Fig. 3 we show the 'errors' in the binding energy using the same 
approximations listed above for the potential. As expected, all error lines drop abruptly to 
zero when A ^ Aq, where the running coupling is chosen to fit what we define to be the exact 
binding energy. Away from this point we can analyze the errors. With the hamiltonian (a) 
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Figure 2: The binding energy for the two-dimensional delta-function potential with various ap- 
proximations for the SRG hamiltonian. The exact theory is fixed by choosing ag = 1.45 at 
Ao = 100. 
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Figure 3: The SRG errors in the binding energy for the two-dimensional delta-function potential 
using various approximations for the SRG hamiltonian. The exact theory is fixed by choosing 
Eo = 1 at Ao = 100. 



(unrenormalized) we obtain a dominant error that scales like In (Aq/A), corresponding to the 
leading order error. With the renornialized hamiltonian (b) the dominant errors scale like 
[ln(Ao/A)]~^, indicating the elimination of the leading order logarithmic errors. With the 
hamiltonian (c) there is a small shift, but no significant change in the error scaling. The 
added irrelevant operator may remove errors of order (Aq/A)^, but these are smaller than the 
remaining [ln(Ao/A)]^^ errors. With hamiltonians (d) and (e) in the range of intermediate 
cutoffs (£"0 << A^ << Ag) there is also only a shift in the errors. The dips in (d) and (e) 
correspond to values of A where the binding energy equals the exact value. 

This behavior is a perturbative artifact that can be understood in the following way. 
Consider the Schrodinger equation with potential (d). Rescaling the momenta p ^ Xp we 
obtain 



P^^iP) - 



and 



Ex = X' 



roc 
JO 



dp p (p^ 



«A,4 

2n 



;^m2 



dq q e 



iP^-Q^) 



Hq) 



^m 



Hp)\ 



"A, 4 
277 



/o~ dpp /o°° dgge(P'-«") l>(p) |.(g) 



/o°° dpp\^ 



(99) 



(100) 



As shown in Fig. 4, the coupling renornialized to fourth-order, 0^,4, approximately freezes 
for small A and as a consequence the bound-state energy scales like Ex — X'^ x constant 
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Figure 4: The SRG running coupling for the two-dimensional delta-function potential obtained 
with a\Q at Aq = 100 fixed to fit -Eq = 1- 



eventually becoming equal to the exact value and then deviating again. With the haniiltonian 
(e) the behavior is similar, with the dip occurring at a different value of A because of the 
irrelevant operator. For small values of A the lines converge, indicating the breakdown of the 
perturbative expansion. 

An alternative prescription is to use the potential derived in subsection 5.1 as the starting 
point for the similarity transformation. We introduce a large momentum cutoff A, define 

0^0=0 = aA = —^ TTv , (101) 



i"(i + f 



and set all irrelevant operators to zero at sq = 0. Note that the coupling a^o is fixed at Aq = 
cx) by fitting the exact binding energy. With this definition the similarity hamiltonian with 
no similarity cutoff becomes well-defined and we can set all of the similarity transformation 
boundary conditions at s = 0. The previous derivation remains essentially the same. The 
only modification is that all integrals over momentum are cut off (p < A). 
At second-order we have 



s,A 



1 da^^^^ 
2tt ds 



.(2) 



1 2 ^^s,A 
— "s.A 



{P,P') 



(2vr) 



ds 



-a 



s,A 



r(2) 



Iliip^p') 



(102) 



where 



tiiip,p') 



{2ny 



dk k (2k'-p^-p'^) e-2.bV^+fc4-fc2{p2+p'= 



20 



1 e^^** f P''' 
(27r)2 4^ 



(103) 



The resulting second-order running coupling and irrelevant operator are given respectively 
by 

as,A,2 



OA 



and 



F^liP.P) 



1 

4 
1 

+ 4 
1 



At third-order we have 



1 - f^ [7 + In (2sA4) - Ei (-2sA4)] 

7 -h ln(2s / p'^) - Ei(-2s / p''^) 
7 + ln(2sA^) - Ei(-2sA^)" 

7 + In (s [(p2 - A2)2 + (p'2 _ a2)2 _ (p2 _ ^.2)2] ) 
_Ei (-5 [(/ _ a2)2 + (p'2 _ ^2)2 _ (p2 _ ^/2^2j ^ 



1 das,A 1 



2ix ds 



(27r)2 "^'^ 



ds 



+ 



(2vr)^ 



3 "s,A 



3 dF^,A^p,p') _ ^2 r(2)/^„^,„,3 .(3) 



ds 



where 



^S(p,p') 



(27r)3 



dk k (2A;2 -p2 _y2) ^-2sbV^+A.4-/c2(p2+p'^)] 



Fi2)(p,A;)+FP(A;,p') 



In this case, 



^i3(0,A:) 



and so 



45(0,0) 



i 
4 


7 + ln(2sA^) 


-Ei(- 


-2sA^)" 




1 
4 


7 + ln(2sA^ - 


- 2sfc2A2) - Ei(- 


2sA^ + 2sfc2A2)] , 


1 

(2vr)3 


1 f^ 

- / dk4k- 

4 7o 


g-2sfc* 


([7 + ln(2sA^)-Ei(-2sA^) 


7 + 


ln(2sA^ - 2sk^A^) - 


Ei(-2sA 


' + 2seA')]) . 



(104) 



(105) 






(107) 



(108) 



(109) 



Since 1^^(0,0) / the term proportional to a^ ^ in Eq. (|85|) does not vanish. To obtain 

(3) \ — I 

as,A,3 we evaluate 1^^(0,0) and solve Eq. (|85| ) numerically. 

(3) 

The equation for F^ jl{p,p') is given by 



1 dF^'l{p,p') 1 



r(2), 



.(2) 



(27r)3 ds 



TT 



iiw,o)F'^'>{p,p')+m{p,p'). 



(3), 



(110) 
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To obtain F^ l^{p,p') the integrals over k and s must be evaluated numerically. In the limit 
sA^ -^ oo with s fixed at some non-zero value 



p{2) 



FILM) 



1 



ii5(o,o)^o, 



2 „/2^ 



7 + ln(2s p^ p^) - Ei(-2s p^ p' 



(111) 
(112) 



and Os^A becomes indeterminate, requiring the coupling to be fixed at some sq ^ 0. In this 
way, we recover the result of the previous prescription. 

Although less trivial, this prescription allows a more transparent error analysis. We can 
also extend the calculation to larger values of the similarity cutoff. A, and analyze the errors 
in a different scaling regime. In Fig. 5 we show the errors in the binding-energy obtained 
using the following approximations for the potential with A = 50: 



(a) marginal operator with coupling (oao), 



Vx,a{p,p) = -—- e P ; 
(b) marginal operator with running coupling renormalized to second-order (aA,A,2) 

VxAp^p') 



I 2 /2x2 



27r 



(113) 



(114) 



(c) marginal operator plus second-order irrelevant operator with running coupling renormal- 
ized to second-order (aA,A,2)-^; 



(2)^ 

A ^' 



VxAp^p) 



aA,A,2 /^"A,A,2 



FikiP^P') 



27r ■ V 27r 
(d) marginal operator with running coupling renormalized to third-order (aA,A,3)) 

VxAp^p') 



ax,A,3 -(zLpl 



2tt 



(115) 



(116) 



(e) marginal operator plus second-order irrelevant operator with running coupling renormal- 
ized to third-order aA,A,3)-^A )' 



VxMp^p') 



aA,A,3 , /Q;a,A,3 



2tt 



+ 



V 2tt 



Fnip^p) 



(117) 



There are clearly two distinct scaling regions when an additional large momentum cutoff 
A is placed on the initial matrix and a similarity cutoff is then applied. When the similarity 
cutoff is larger than A, we see power-law improvement resulting from the addition of irrelevant 
operators. Curves (a), (b), and (d) all have the same slope. None of these hamiltonians 
contains irrelevant operators, but the marginal coupling differs in each. All results become 
exact as the similarity cutoff goes to infinity, and these curves are close to one another because 
the coupling runs little in this region. Curves (c) and (e) show that there is a power-law 
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Figure 5: The SRG errors in the binding energy for the two-dimensional delta-function potential 
using various approximations for the similarity hamiltonian. The exact theory is fixed by regu- 
lating the "bare hamiltonian" using a sharp momentum cutoff, A, and letting the bare coupling 
depend on A such that the binding energy is fixed. We use A = 50 and ii^o = 1- 
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improvement when irrelevant operators are added, and that once again when the similarity 
cutoff is larger than A an improvement in the running coupling makes little difference. Even 
though the coupling in front of this operator is approximated by the first term in an expansion 
in powers of the running coupling, the coupling is sufficiently small that this approximation 
works well and the operator eliminates most of the leading power-law error in curves (a), 
(b), and (d). 

When the similarity cutoff become smaller than A we see a crossover to a more complicated 
scaling regime that resembles the SRG scaling discussed above. The error displayed by curve 
(a) approaches 100%, while the running coupling introduced in curve (b) reduces the error 
to an inverse logarithm. Improving the running coupling in curve (d) further reduces the 
error, and we see that curve (c) crosses curve (d) at a point where improving the running 
coupling becomes more important than adding irrelevant operators. As above, the best 
results require us to both improve the running coupling by adding third-order corrections and 
add the second-order irrelevant operators. In no case do we achieve power-law improvement, 
because as we have discussed there are always residual inverse logarithmic errors. Had we 
fit the running coupling to data, as we would do in a realistic calculation, we would obtain 
power-law improvement and the residual error would be proportional to an inverse power of 
the cutoff times an inverse power of the logarithm of the cutoff. 

In Fig. 6 we show the running coupling at 2nd and 3rd order. Although the 3rd order 
corrections are small for all A and vanish when A ^ oo, the improvement resulting from this 
correction in Fig. 5 is significant. 

The scaling behavior with a large momentum cutoff A in place is complicated, but it 
is fairly straightforward to understand it and to find a sequence of approximations that 
systematically improve the non-perturbative results. The calculations become increasingly 
complicated, but at each order one must improve the running coupling, or fit it to data, and 
add higher order irrelevant operators. In a field theory we need to let A ^ cx) and study 
the scaling behavior of the theory in the regime where A ^ A. Although we do not display 
a compete set of figures, in Fig. 7 we show what happens to the running coupling as A is 
increased, with the bound state energy fixed at one. 

As is evident in the exact solution, as A increases the coupling decreases. When A ^ A, 
the coupling runs slowly and stays near its asymptotic value. As A approaches A the coupling 
begins to run noticeably, and when A becomes much less than A the coupling approaches a 
universal curve that is insensitive to its asymptotic value. Plots of the error in the binding 
energy for various approximations and different values of A closely resemble Fig. 11, with 
two scaling regimes whose boundary is A = A. 

We close this section by reminding the reader that in all of these calculations there is only 
one free parameter. In a realistic calculation we would fit this parameter to a binding energy 
and we would expect to see residual errors in other observables that is inversely proportional 
to powers of the cutoff and logarithms of the cutoff. 

4 Conclusions 

We have illustrated the similarity renormalization group method for producing effective cutoff 
hamiltonians using the two-dimensional delta-function potential. We have shown that the 
SRG with coupling coherence leads to errors that scale as inverse powers of the cutoff and 
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Figure 6: The SRG running coupling for the two-dimensional delta-function potential renormalized 
to second and third-order obtained with aAo=oo -^ ola = 47r/ln fl + ^j. 
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Figure 7: The SRG running coupling for the two-dimensional delta-function potential renormalized 
to second-order obtained with aAo=oo — * c^a = 47r/ln f 1 + ^j. 
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inverse logarithms of the cutoff. The SRG with couphng coherence requires the same number 
of parameters as the underlying 'fundamental' theory, but the cost is exponentially increasing 
algebraic complexity to remove errors that contain inverse powers of logarithms of the cutoff. 
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